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Abstract 

We propose a new version of the ETAS model, which we also analyze theoreti¬ 
cally. As for the standard ETAS model, we assume the Gutenberg-Richter law 
as a probability density function for background events’ magnitude. Instead, the 
magnitude of triggered shocks is assumed to be probabilistically dependent on 
the triggering events’ magnitude. To this aim, we propose a suitable probability 
density function. This function is such that, averaging on all triggering events’ 
magnitudes, we obtain again the Gutenberg-Richter law. This ensures the va¬ 
lidity of this law at any event’s generation when ignoring past seismicity. The 
probabilistic dependence between the magnitude of triggered events’ and the one 
of the corresponding triggering shock is motivated by some results of a statistical 
analysis of some Italian catalogues. The proposed model has been also theoreti¬ 
cally analyzed here. In particular, we focus on the inter-event time which plays 
a very important role in the assessment of seismic hazard |7]. Using the tool of 
the probability generating function and the Palm theory, we derive the density 
of inter-event time for small values. 


1 Introduction 

The Epidemic Type Aftershock Sequences (ETAS) is a model based on a specihc 
branching process and represents a benchmark in statistical and mathematical seis¬ 
mology. The process generates the events along consecutive generations. At any step, 
each event of a given generation produces, independently of the others, its own cluster 
of events. The number of elements in a cluster is a Poisson random variable dependent 
on the generating event’s magnitude through the productivity law. Occurring times 
and magnitudes of each event in a cluster are independent of each other and are dis¬ 
tributed according to the Omori-Utsu and the Gutenberg-Richter laws, respectively 
(for the above-mentioned three laws see equations fl72]) . flT^ and fl71]) in Appendix lA)) . 
These times and magnitudes are also independent of the characteristics of past events. 
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Similarly, given an event, the elements relative to its cluster are spatially distributed 
according to a particular model centered at it (see 0 )- Here we will consider the pro¬ 
cess only regarding to the time, ignoring the spatial location of events. In particular, 
we will focus on the random variables relative to inter-event times, that are the times 
between consecutive shocks. 

In the ETAS model, seismic events are divided into spontaneous (or background) 
and triggered ones. The characteristics of the branching process involved and the ho¬ 
mogeneity of the background component imply that the total process is homogeneous, 
too. Therefore, we can deduce that inter-event times have all the same distribution. 
We will then consider the single inter-event time between consecutive shocks. The 
characteristics of this random variable play a very important role in the assessment 
of seismic hazard [ 3 . In what follows, we will assume that A is the rate of the whole 
process of’’observable” events and u is the one corresponding to the ’’observable” events 
of the spontaneous component. 

Observation 1. Given the events of any seismic catalogue, each of them is consid¬ 
ered observable only if it has a magnitude larger than a certain cutoff value m. This 
threshold is known as the completeness magnitude^ that is the smallest value such 
that catalog events’ magnitudes exceeding that value are distributed according to the 
Gutenberg-Richter model. The deviations from that model, when considering cata¬ 
logue magnitudes smaller than m, can be interpreted in terms of wrong or missing 
events measurements. 

A very useful tool to study the distribution of the inter-event time is the probability 
generating function (PGF) of the random number N{t) of observable events occurring 
in the time interval [0, r], that is the quantity 

GMir){z)=nZ^^^\ ( 1 ) 

The value of this function in zero gives the probability P(t) to have zero events in 
[0,r]. The inter-event time density is then obtained through the Palm equation by a 
double time derivative of P(t) and scaling. This has been done for the standard ETAS 
model by Saichev and Sornette [T3]. However, the standard ETAS model assumes that 
the magnitude of each event in any cluster is distributed according to the Gutenberg- 
Richter law, independently of any past event. Instead, it is assumed that the number 
of triggered events in a cluster depends on the triggering event’s magnitude. We hnd 
intuitive and more realistic to also assume that the distribution of triggered events’ 
magnitude in a cluster depends on the triggering events’ magnitude. This is supported 
by evidences of an empirical study (see Cl)- Therefore we propose a new ETAS 
version which assumes the above-cited hypothesis. Following the same technical tool 
of Saichev and Sornette, that is the above-mentioned PGF, we derive an expression 
for the inter-event time density for small values. In the following section we are going 
to compute the above-mentioned probability generating function corresponding to the 
proposed model. 
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2 Probability generating function for earthquake se¬ 
quences 

As previously specified, we are interested in N{t), that is the number of observable 
events in [0,r], which are of two types: spontaneous and triggered. Let’s consider at 
first the triggered events. We can associate to each of these events the spontaneous 
shock from which the branching process reaches that event. Let’s consider then all 
spontaneous events. The triggered events can be partitioned according to their asso¬ 
ciated spontaneous triggering shocks. The total number Ni{t) of observable events 
triggered in [0, r] will be the sum of observable events triggered by each spontaneous 
event. Let’s consider now a partition of the infinite interval of times smaller than r 
into small subintervals Ik = [tk,tk + At] grouped in two categories. To the first (sec¬ 
ond) category belong subintervals of kind Ik^ C (—oo,0) C [0,r]). The probability 
generating function of the number Ni{t) can be obtained as the limit, with respect 
to At, of the corresponding probability generating function relative to the subinterval 
partition. 

Recall now both the property of the PGF expressed in equation (I76p . Appendix [Bl 
and the fact that there’s independence between the numbers of events, triggered by the 
spontaneous ones belonging to disjoint subintervals. We can therefore deduce that the 
PGF of Ni{t) is the limit of the product between the respective probability generating 
functions relative to the different subintervals. Let’s consider the PGF of the number 
of events in [0,r], triggered by spontaneous shocks contained in a generic subinterval. 
Thanks to the properties of the probability generating function, this PGF is obtained by 
composing the PGF of the number of background events occurred in that subinterval, 
with the PGF relative to the number of events in [0, r] triggered by just one of the 
considered spontaneous events (see Proposition fl75|) . Appendix iBj) . 

Furthermore, in order to obtain the probability generating function of the total 
number N{t) of observable events, we have to take into account spontaneous events, 
too. Obviously, the background events to be considered are only those which fall in 
subintervals of the second type. In this case, for each of these events we will count 
all the shocks triggered by it, plus the event itself. This last is equivalent to multiply 
by z the probability generating function of the number of events triggered in [0, r] by 
a single spontaneous event of magnitude m', itself belonging to [0,r]. More precisely, 
since this event has to be observable, the multiplication by z will be done only if its 
magnitude is larger than the threshold m. Therefore, instead of the multiplication by 
2 :, we will multiply by 


f{z,m\m) = U{m'<m)-IzU{m'>rn), (2) 

where U{x > 0) is the indicator function having value one where x > 0, zero otherwise. 

Let’s now compute explicitly the probability generating function G(-;t) of N{t). 
Let’s consider the number of events with magnitude larger than m, triggered in a 
generic interval [0, zz] by a background event with magnitude m' occurred in a time 
—s, where s > 0. In what follows, we will indicate with G{--, s,u,m') the PGF of 
such a number. Recalling that the random variable associated with the number of 
spontaneous events in a generic subinterval of length At follows a Poisson distribution 
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with parameter uAt, for the previous reasonings we have 
r) = 


lim ]^exp< —uAt[l — G{z; s,T,m')] > 

4i I J 

• exp<|— C(;At[l — f{z, m!, m)G{z; 0,t — w, m')] |> 

[1 — G{z] —t',T,m')]dt' 


= exp < —u 


+ / [1 — f{z,m',m)G{z-,0,T — t',m')]dt' 


= exp < —u 


[1 — G{z-, t, r, m')]dt 


+ J [1 — f{z,m',rn)G{z]0,t,m')]dt 


( 3 ) 


where — s G and w ^ Ik^- 

The formulae obtained till here are relative to observable events with magnitude 
larger than m, triggered by a spontaneous shock with a hxed but generic magnitude 
equal to m'. Nevertheless, we have to introduce another threshold value mo, called 
reference magnitude, which corresponds to the minimum value for an event to be able 
to trigger their own offsprings [6]. We therefore have to substitute in ([3]) the generating 
function G{z] •, •, m') with 



■,m')p{m')dm', 


where p{m') is the Gutenberg-Richter law (equation (1721) in Appendix Rl). 

Observation 2. Reference magnitude is usually set less than or equal to the complete¬ 
ness one ng, Indeed, even if the opposite inequality is verihed, we could however 
reduce to the previous case by excluding the data relative to observable events with 
magnitude less than mo- 

The function p{-) is universally adopted to model the probability density function 
of a generic event’s magnitude, if we don’t take into account the characteristics of 
its previous seismicity. In the specihc case we are considering background events; we 
assume that their magnitude is also modeled by the function p{-). Obviously it holds 


Setting 


p{m')dm' 


'mo 


1 . 


L{z;t) 



[1 — G{z] t, r, m')]dt 


+ J [1 — f{z,rn',m)G{z;0,t,m')]dt^ , 


( 4 ) 
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we can write, with abuse of notation, 


n{z-T) = 


( 5 ) 


Let’s substitute now the above-cited function f{z,m',m) (equation ([2])) in the last 
expression of (jl]). We obtain 


' mo 


L{z;t) = / dm'p{m')i / [1 — G{z-,t,T,m')]dt 


1 — [U(m' < m) + zU{m' > m)]G{z-, 0, t, m!) 


dt 


'0 


G{z] t, T, m')p{m')dm! 


'mo 


dt 


1— / G{z]0,t,m')p{m')dm' 

/ 0 L J mo 

poo 

— z G{z-,0,t,m')p{m')dm' 

J m 


dt 


1 — / G{z] t, T, m')p{m')dm' 


'mo 


dt 


G{z; 0, t, m!)p{m!)dm! 


'mo 

^OO 


-I-(1 — z) / G{z;0,t,m')p{m')dm' 


dt. 


( 6 ) 


At this point, we want to compute explicitly the probability generating function 
r2(-; r). To do that, we hnd an expression for G{z-, t, r, m') that is the PGF of the num¬ 
ber of observable events in [0, r] triggered by a hxed spontaneous shock with magnitude 
m' occurred in —t, where f > 0. Let’s consider now any observable event triggered in 
[0,r] belonging to any generation from the second one on and let N 2 {t) denote their 
total number. As before, we can associate it with a hrst generation triggered shock 
occurred between the spontaneous one and the event itself. The number N 2 {t) will be 
the sum of observable events triggered by each hrst generation triggered shock. Thanks 
to the auto-similarity of the process, the above-mentioned probability generating func¬ 
tion G{z] •,•,•) is the same for both a background event and a triggered one. Therefore, 
as previously done, we obtain 


G{z-,t,T,m') = exp<| <h(—s-|-f)Afp(m')[G'(z; s, r, m") — 1] 


4^c[-i,o] 


JJ^expJ $(t(; -I- t)Atp{m')[f{z, m", m)G{z] 0,t — w, m") “ 1] p (V 


lk2 


where —s E Iki, w E and m" is the magnitude of the hrst generation triggered shock. 
Recall that <h(-) is the Omori-Utsu law and p(-) is the productivity one (respectively 
equations fl73|) and frM|) in Appendix]^. Again, we will have to average with respect to 
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the first generation triggered shock’s magnitnde. To this end, we assnnie that the latter 
magnitnde, here indicated with m" conditioned on the backgronnd event’s magnitnde 
m', follows the model P(m"|m') = p{m''\m'), where the proposal for the last probability 
density is given in Appendix O Obviously it must hold 


p{m"\m')dm" = 1. 


' rno 


Any hrst generation triggered event considered is also a triggering one. Then, its 
magnitude will have to be integrated from mo to inhnity. Recalling equation (]2|), 


6 


from (jT]) we obtain, with abuse of notation, the following expression: 


POD 


G(z; f, r, m') = exp< —p(m') / dm"p{m"\m') 


' mo 


<h(f + x)[l 


’-t 


G( 2 :; —x, r, + / <h(f + x)[l 


f{z, m”, m)G{z] 0, r — x, m'')\dx 


= exp< —p(m') / dm"p{m"\m') / <h(f + x)(ix 

I J mo i 

/ O nr 

<h(f + x)G'( 2 ; —X, r, m")(ix + / <h(f + x)(ix 
-i Jo 


f{z,m",m) / + x)G(z; 0, r — X, m")(ix 


= exp< —p{m!) / dm''p{m''\m!) 


' mo 


$(x)(ix 


LJo 


ft ft+T 

<h(f — x)G'(z; X, r, m")(ix + / <h(x)(ix 


'0 


f{z,m",m) / + r — x)G(z; 0, X, m")(ix 


/o 

poo 


= exp< —p(m') / + r) 

[ Jmo 

- (^'(•) 8)G'(^;-,r, m"))(t) 

- f{z, m", m)(<h(f + •) O 0, •, m"))(r)] 

= exp<[—p(m') 6(f + r) 

($(•) 0 G(x; •, r, m"))(f)p(m"|m')(im" 

) 

($(f + •) 0 (^(z; 0, •, m"))(r)p(m"|m')(im" 

poo 

+ (1 — 2 :) / (<h(f + •) 0 G( 2 ;; 0, •, m"))(r)p(m"|m')(im" 


' mo 

POO 


' mo 


where 0 is the convolution operator and 


b{t) = / ^{x)dx. 

Jo 


( 8 ) 

( 9 ) 
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At this point, in order to find the probability of having zero events in [0, r], we will 
have to snbstitute eqnation ([8]) in ([6]) opportunely, reminding that 

P(r) =fi(0;r) = 

In what follows we will work separately on two cases. First of all on the particular case 
m = rriQ, then on the general one m > mo. 


3 The case m = mo 

In this section, we want to focus on the particular case in which m = mo- We will 
focus on this case because we can find a closed-form expression for the probability 
P(r) of having zero events in [0,r], for small values of r. Than, by using the Palm 
theory, we derive a closed-form expression for the density of inter-event time. The 
relevance of the case m = mo can be understood by the following reasoning. On 
one side, the reference magnitude mo is connected to the intrinsic physical features 
of soil and remains almost constant in a sufficiently long time frame. On the other 
side, the completeness magnitude is related to the number of seismic stations and their 
sensitivity. Then, the latter magnitude varies with seismic stations density. Hoping 
for a gradual improvement of technology both in instrumental sensitivity and in the 
number of seismic stations available, we can certainly suppose to reduce m till it will 
be just equal to mo- 

Setting m = mo in the last expression of formula ([2]) and recalling (E]), we find that 
the probability generating function r) relative to the total number of observable 
events in [0,r] is 


17 ( 2 ;; r) = exp< —uj 


1 — / G{z; t, T, m')p{m')dm' 


'mo 


dt 


+ 


'0 


1 — z G{z-,0,t,m')p{m')dm‘ 


'mo 


dt 


In this particular case the function L(0; r) becomes 


L(0;r) = 


'mo 

roo 


' mo 


G(0; t, r, m')p{m')dm' 
G(0; t, r, m')p{m')dm' 


dt+ dt 


dt + T 


Nm=moit,r)dt + T, 


with 


Nm=mo {t,T) = 1- I G(0; t, r, m')p{m')dm' 


' mo 


( 10 ) 


( 11 ) 


( 12 ) 


Let’s compute then the probability generating function G(0; t, r, m'). Setting z = 0 
and m = mo in the last expression of ([S]) and considering the approximation of the 
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exponential function, we obtain 


G{0;t,T,m') =exp< —p{m') 


h{t + r) 


($(•) (8) G(0; •, r, m'')){t)p{m''\m')dm'' 


' mo 


— p{m') 


b(t + t) 


($(•) 0 G(0; •, r, m'')){t)p{m''\m')drri'' 


'mo 


(13) 


This approximation in Taylor series around the point zero makes sense since the expo¬ 
nential argument is zero when r = 0. In fact, G(0; •,0,m") = 1. If we substitute the 
previous formula in (IT^ . we get 


m=mo(ty 


I mo 

f*00 


dm'pim') 1 1 “ pijn') b(t -|- r) 

— I ($(•) 0 G(0; •, r, 

J rtiQ 

poo 

=1 — 1 + b{t + t) / p{m')p{m!)dm' 

J mo 


dm'p{m')p{m!) / ($(•) 0 G(0; •, r, 


'mo 


'mo 


=nb{t + t) — / dm'p{m')p{m') / dm"p{m"\m') 


' mo 


' mo 


■ / ^{x)G{0;t — x,T,m")dx 

Jo 

switch the order of integration 

pt poo 

=nb(t + t) — / dx^{x) / dm''G{Q]t — x,T,m") 
Jo J mo 


p{m!) p{m')p{m'' \ m')dm' 


' mo 


=nb{t + r) — / dx^{x) / G{0;t — x,T,m")lA{'m")dm'', (14) 

Jo J mo 


with 


lA{'rn") = / p{m')p{m')p{m''\m')drn! 

J mo 


(15) 
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and 


n = 


p{m')p{m')dm' 

poo 

= /3k 


'mo 

roo 


'mo 


'mo 


13k 


-(/5 - a) 

13k 


^-(P-a){m'-mo) 


mo 


(3 — a 


(16) 


Equation (fT6|l has been obtained using the formulae fl72|) and fl7^ in Appendix 
Let’s assume now that the following condition is verified: 


p[m') p{m')p{m''\m')drn' 


p{m')p{rn')dm' 


= p{m" 

t 


I Aim") = np{m"). 


Observation 3. This condition represents the invariance of the Gutenberg-Richter law, 
weighted by p{m'), with respect to the conditional probability density function as¬ 
sumed valid for triggered events’ magnitude. It is a focal point of our study. Although 
we assume the existence of a density of triggered events’ magnitude different from 
the Gutenberg-Richter law, at the same time the above condition justifies its validity. 
More precisely, the condition we are imposing tells us that, if we average the conditional 
density p{m"\m') over all the possible triggering events’ magnitudes m', distributed ac¬ 
cording to the Gutenberg-Richter law and taking into account the productivity model, 
we obtain again the above-cited law. This shows that first generation events’ magni¬ 
tude follows the Gutenberg-Richter law as assumed for background shocks. In addition, 
this property will be obviously true for events’ magnitude of any generation. Hence, 
iterating the reasoning, a generic event’s magnitude, whatever generation it belongs 
to, will be distributed according to the Gutenberg-Richter law. 
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Under this assumption, equation ffl^ becomes 


Nm=mo it) 't) ^nb{t + t) — n I dx <h(a:) / G(0; t — x,t, m")p{m")dm" 


' mo 


=n[b{t + r) — b{t)] + n 


b{t) 


dx^{x) I G{0-,t — x,T,m")p{m")dm" 


' mo 


=n\b(t + t) — b(t)\ + n / dx^{x) 


^(0; t — x,T, m")p{m")dm" 


f mo 


=n\b{t + t) — b{t)~\ + n / ^{x)Nm=mo{t ~ ^)T)dx, 


(17) 


where in the last equality we’ve used formula 0121) . For simplicity of notation, in what 
follows we will substitute approximation symbol with equal sign. 

In order to obtain the expression to substitute in flTT]) . let’s compute the temporal 
integral of both the hrst and the last members of equation flTT)) . We get 


r’OO pt 


Nm=mo{t)T)dt=n j [b{t + t) - b{t)\dt + n j 

Jo Jo Jo 

switch the last two temporal integrals 


- x)Nm=mo (x, T)dxdt 


/*oo 


=n / [b{t + r) — b{t)~\dt + n / / — x)Nm=mo{x,T)dtdx 

Jo Jo Jx 

POO POO 

=n / [b{t + t) — b{t)~\dt + n / Nm=mo{x,T)dx, 

Jo Jo 

where the last equality is obtained because, recalling 075)) . it obviously holds 

POO POO 

/ — x)dt = / ^{y)dy = l. 


( 18 ) 


From flTSj) we then have 


Nm,=mo (t) X)dt 


n 


1 — n 
n 


[b{t + r) — dt 
^0 


'0 


1 — n 

nc^ r 


1 - 


{c + t + tY 
1 1 


-1 + 


^ ^ Jo [{c + ty {c + t + ty 


(c + ty 
dt 


dt 


nc 


{c + ty ^ - (c + t + r) 


{l-n){l-e)l 

n 

nc 

{l-n){l-e) L 


i-e 


(c + t) 


i-e _ ^1-9 


(19) 
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The second equality follows from the explicit form of the function b{t) dehned in dOj): 


b{t) = 


9c^ 


(c + 
rt 


dx 


= 9r (c + X 


,-1-0 


= -c^(c + x)-^\i 


= 1 - 


(c + ty 


On the other hand, the last equality in (IT^ is obtained since the limit of the function 
(c + ty~^ — (c + f + is, for t tending to inhnity, always zero. In fact, for ^ > 1 it 
is obvious; if instead 1 > 0 > 0 we have 


lim r(c + ^ — {c + t + tY = lim (c + t) 

L -I t^OO 


1-0 


t^OO 


1-1 + 


r 


c + t 


i-e' 


1-1 + 


= lim 


c+t 


1-0 


t^oo (c + ty ^ 
apply De L’Hopital’s rule 




c+q2 


t^oo (0 — l)(c + f)^ ^ 
-e 


t^oo (c + ty 

= 0 . 


At this point we can substitute equation flT^ in flTTD : 


L(0;r) = 


Introducing the functions 


n 


(l-n)(l-0)L 


y{c + r)^ ^ — c 


+ r. 


a{t) = 1 — b{t) = / ^{x)dx = 


and 


we get 


— / a{x)dx = - -x(c + t) 

./n 1 — 9 


\i-0r 
lo 


{c + ty 

[c®(r + cy~^ — c] 


frl = 


= exp< —UT 


noj 


1 — n 


A(t) 


( 20 ) 


( 21 ) 


( 22 ) 
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As already specified, using the Palm theory we obtain the density iA('r) relative to the 
inter-event time in the following way (see ID): 


1 d2p(r) 

A 


(23) 


Let’s compute than the second derivative of the above-mentioned probability. We have 


in fact 


Hence, 


dP(T) 

dr 


= P(r) 


nu 


-u 


n 


■a(r) 


d . d r f f . 

—A[t) = — a[x)dx = a[T). 
(IT dr Jq 


(24) 


dr'^ 


=P(r) 

=P(t) 


-u 


00 + 


nuj 
1 — n 
noj 
1 — n 

noj 


a(r) +P(r)\^{c + r) 


\-e-i 


1 — n 


= exp < —UJT — 


1 — n 


ar + 


A{t) 


noj 
1 — n 


$(r) 


00 + 


nu 
1 — n 


a r -h 


noo 
1 — n 


<h(r) 


In conclusion. 


^ / N 1 I "noo . 

Frir) =Yexp<^ -ut - - - A{t) 

A 1 — n 


noo / noo . 

00 + - -afr) ) -F - $(t) 


1 — n 


1 — n 


(25) 


We can see that this density doesn’t depend on the assumed model for the conditional 
distribution of triggered events’ magnitude. Therefore, in the case m = mo, in which 
all the events are observable, our hypothesis doesn’t play any role. 

4 The general case m > mo 

In order to hnd the expression of the inter-event time density for the general case, let’s 
start again with formula (E]). Let’s define the function 


>I'(;i('), in) = / f)(m/)e-'’*'"'>'*<'"')*„'_ 


(26) 


where m > 0 and h{-) is a continuous function dependent on some variables; it is 
here considered as a function of only one of them, corresponding to the triggering 
event’s magnitude m'. Recalling expression (|8]) of the probability generating function 
G{z; t, r, m'), we can rewrite the integral G{z; t, r, m')p{m')dm' of the last member 
of (El) as 


G{z;t,T,m')p{m')dm' = '^[y{z;t,T,mo,m'),mo], 


(27) 


' mo 
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where 


y{z-,t,T,mo,m') = b{t + t) — / dx^{x) / G{z;t — x,T,m'')p{m''\m')dm'' 


'mo 


dx^(t + r — x) I G{z;0,x,m'')p{m"\m')dm'' 


'mo 


+ (1 — z) / dx^{t + T — x) / G{z]0,x,m'')p{m''\m')dm‘ 

Jo J m 

=b{t + r) - ($(•) ® D{z; •, r, mo, m'))(t) 

- (<h(t + •) (8)-D+(z; •,mo,m'))(r) 


+ (1 — 2;)($(t + •) 0 D^i^z] 

•,m,m'))(r). 

(28) 

In the latter equality we have set 



poo 

D{z;t,T,mo,m) = / G{z;t 

, r, m')p{m'\m)dm', 

(29) 

Jmo 



poo 

D^{z;t,m,m) = / G{z;0, 

t, m')p{m'\rn)dm'. 

(30) 


Let’s notice that in our case we have m = rriQ, m. 

Similarly, the other two integrals with respect to magnitude 


G{z;0,t,m')p{m')dm' and / G{z;0,t,m')p{m')dm', 


I mo 


in the last member of ([6]), can be expressed as 

roo 

/ G{z] 0, t, m^)p{m^)dm' = t, mo, m'), mo], 

Jmo 

poo 

/ G{z;0,t,m')p{m')dm' = t, mo, m'), m], 

J m 

where 


(31) 


s{z;t,mQ,m') =b{t) — / dx^{t — x) / G{z;0,x,m'')p{m''\m')dm" 


'mo 


+ (1 — 2 ;) / dx^{t — x) / G(z; 0, X, m )p(m |m')(im' 
Jo J m 

=b{t) - ($(•) 0 D+{z-, •,mo,m'))(t) 

+ (1 - 2 ;)(<h(-) 0 Dj^{z-, •,m, m'))(t). 


(32) 


Concluding, equation ([H]) becomes 

L{z]t) = j |l - ^[y{z;t,T,mo,m'),mo\^dt 

+ f |l — \h[s( 2 ;; t, mo, m'), mo] + (1 — 2 ;)'h[s( 2 :; t, mo, m'), m]|(it. (33) 
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As explained before, in order to find the probability P(t) of having zero events in a 
temporal interval of length r, we have to compute the probability generating function 
in zero. We then have 

P(r) = (34) 

Obviously, the expression of the function L(0; r) can be obtained by setting z = 0 
in fl33|) . In particular we get 


L(0;r) 



1 — 'h[s(0; t, mo, m'), mo] + 'h[s(0; f, mo, m'),m] >dt 


N{t,T)dt+ / N_{t)dt, 


(35) 


where 

N{t, t) =1 - 4^[?/(0; t, T, mo, m'), mo] (36) 

and 

iV_(f) =1 — d^[s(0; t, mo, m'), mo] + 'h[s(0; t, mo, m'), m]. (37) 

Observation 4. Recalling equations fl27|) and fIM]) . it is clear that we can rewrite in the 
following way the last two functions just dehned: 

poo 

iV(f,r) = l— / G{0;t,T,m')p{m')dm' (38) 

J mo 

and 

poo poo 

N_{t) = l— / (7(0; 0, f, m')p(m')dm'+ / (7(0; 0, f, m')p(m')(im'. (39) 

J mo ^ m 

Obviously, in the case m = mo we have that N(t, r) = Nm=mo{t, 't) and N_{t) = 1. 

At this point, in order to obtain the probability P(t) of having zero events in 
[0,r], we have to hnd an expression for the functions iV(t, r) and N-{t) and then 
substitute the obtained results in equation fl35|) . Actually, since they are integrals 
rather difficult to compute, we can approximate them by means of the hrst order Taylor 
power series expansion of the function T(/i(-),m), with respect to the hrst argument, 
around the point zero. In our case it holds m = mo,m and, recalling equations fl28|) 
and (l32|) . the function h{-) in iV(f, r) and iV_(f) is respectively equal to y{0;t,T,mo, •) 
and s(0; t, mo, ■), where 

?/(0; t, T, mo, m') =b{t + r) - ($(•) (g) i7(0; •, r, mo, 

- {^{t + •) ® T)+(0; •, mo, m'))(r) 

+ (<I)(t + •) ® D+(0; •,m,m'))(r) (40) 
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and 


s(0;t, mo,m') =b{t) - ($(•) (g) I^+(0; •,mo,m'))(t) 

+ ($(•) (g) D+(0; (41) 

We notice that the previous approximation makes sense since, in both cases, the hrst 
argument of 4/ is small when r ~ 0. In fact, we obtain that 

• for the argument holds 

1 /( 0 ; f, r = 0, mo, m') =b{t) — b{t) = 0. 


This follows because, recalling that G{z-^ t, r, m") is the probability generating 
function of the number of events triggered in [0,r] by a shock with magnitude 
m', occurred t seconds before time 0, we have 


D(0; f, r = 0, mo, m') = / G(0; t, r = 0, m")p(m"|m')(im" 

J mo 

poo 

= / [0°P(A^(0) = 0) 

^mo 

+ 0^P(A^(0) = l)]p(m"|m')dm" 


The previous calculus is valid because of the regularity of the process and the fact 
that the probability P(A^(0) = 0) of having zero events in an interval of length 
zero is equal to one. 

• on the other hand, for s(0; f, mo, m') holds 

s(0;t = 0,mo,m') =0, 
because t varies between 0 and r = 0. 


Let’s compute then the above-cited expansion and the relative approximations of 
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N{t,T) and Let’s start with the hrst function: 

N{t, t) =1 - ^[y{0] t, T, mo, m'),mo] 

=1 - 


'mo 


' mo 

POO 


' mo 


POO 

- / p{m')p{m')y{0;t,T,mo,m')dm' 

Jmo 

pijri)p{m' )y{Q-, t, r, mo, m')dm' 
p{m')p{m')dm'b{t + r) 

•*oo 

p(m')p(m')(<h(-) 0 D(0; •, r, mo, m')){t)dm' 
p(m')p(m')(<h(t + •) 0 -D+(0; •, mo, m'))(r)(im' 


' mo 

/•CXD 


'mo 

POO 


+ / p(m')p(m')(<h(t + •) 0 L)+(0; •, m, m'))(r)(im', 

J mo 


(42) 


where we have used equations fl26|) . fl36|) and (HUj) . We are going to study separately 
the integrals in the last equality of fH2|) . In the matter of the hrst, recalling the hrst 
equality in f[T6|) . by dehnition we have 


p{m')p{m')dm'b{t + r) = nb{t + r). 


'mo 


For the second one, recalling fllSp and fl2^ . we have 

POO 

/ p(m')p(m')(<!)(•) 0 D(0; •, r, mo, m')){t)dm' 


'mo 


dm!p{m!)p{m!) / rfx $(t — x) G(0; x, r, 


' mo 
rt 


' mo 


f*00 POO 


= / dx^{t — x) 


p{m')p{m')G{0-, X, r, m")p{m"\ m')dm"dm' 


'mo Jmo 

POO 


f*00 


= / dx^{t — x) / dm"G{0]X^T^m") / 


'mo 

POO 


mo 


= / —a:) / G{^-,x,T,m")lA{'m")dm". 


'mo 


If we furthermore use equation ([30]), we can compute the third integral of the last 
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expression in (H2|) as follows: 


p{m')p{m'){^{t + •) ® -^+(0; *, mo, m'))(T)(im' 


'mo 


dm'p{m')p{m^) I dx^{t + T — x) I G(0; 0, x, m")p(m"|m')(im" 


'mo 


'mo 


f*00 /‘CO 


= / dx^{t + T — x) 


p{m')p{m')G{0; 0, x, m'')p{m''\m')dm''dm' 


'mo ^mo 
poo 


= / (ix$(t + r —x) / dm"G{0;0, x,m") / p{m')p{m')p{m''\m')dm' 


' mo 
poo 


mo 

//\ T ( _ ft\ J _// 


= / (ix$(t + r —x) / G( 0 ; 0 , X, m )/yi(m )rfm' 

^0 J mo 

Similarly, the fourth integral becomes 

POO 

/ + •) ® -D+(0; •, m, m')){T)dm' 


'mo 


dm'p{m')p{m') / dx^{t + r — x) / G(0; 0, x, m")p(m"|m')(im" 


'mo 


poo poo 


= / dx^{t + T — x) 


p{m')p{m!)G{Q] 0, X, m")p{m''\m')dm"dm' 


' mo ^ m 
poo 


= / (ix$(t + r —x) / 0, X, m") / 


' m 
poo 


mo 

= I dx^{t + r — x) I G{0-,0, x,m")lA{m")dm". 


We can then hnd the approximated expression of the function N{t,T) by substituting 
in equation fl^ the obtained results. We get 


N{t,T) ^nb(t + r) — I dx^(t — x) I G{0-, x,T,m")lA{m")dm" 


'mo 


dx^{t + T — x) I G{0;0,x,m")lA{'>Ti")dm" 


’ mo 
poo 


+ / dx$(f + r —x) / G{0;0,x,m”)lA{Tn")dm" 


(43) 


For simplicity of notation, from now on we will substitute approximation symbol with 
equal sign. 

Before proceeding with the study of the approximation of the second function N_ (t) 
in fl35|) . we focus for a moment on the approximation obtained for iV(f,r). As we’re 
going to see, it can be demonstrated that 


n 


N{t, T)dt = - -(a(-) O A^_(-))(r). 

1 — n 


(44) 
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Consequently, the function L(0; ■) in fl35|) can be written in terms of the function A^_( 
only. We have indeed 


L(0;r)= / N{t,T)dt+ / N_{t)dt 


n 


1 — n 


:a(-)®iV_(.))(r) 


N_{t)dt. 


(45) 


In order to demonstrate equation fl44|) . we impose in fHHl) that is lA{m") = np{m"), 
as already assumed in our model. This condition, as explained in Observation [3l 
guarantees the validity of the Gutenberg-Richter law when averaging on all triggering 
events’ magnitudes. We then get 

pt poo 

N{t,T) =nb{t + t) — n / dx^{x) / G{0;t — x,T,m")p{m")dm" 


'mo 

^oo 


n I dx + T — x) I G{0-,0,x,m")p{m")dm" 


' mo 
roo 


+ n dx^{t + T — x) / G{0]0,x,m )p{m )dm‘ 


=n[b{t + t) — b{t)] + n< / dx^{x) 


G(0; t — x,T, m")p{m")dm' 


' mo 


+ n< / dx^{t + T — x) 


G(0; 0, x, m")p{m")dm‘ 


' mo 


+ 


G(0; 0, X, m )p{m )dm: 


f m 


=n / ^{y)dy + n / ^{x)N{t — x,T)dx 


n 


rt+r 


dy^iy) 


G{0; 0,t + T — y, m")p{m")dm" 


' mo 


+ / G(0; 0, t + r — 7/, m )p(m jdm' 

J m 

/ t^T 

dy^{y) 

poo 

— I G{0-,0,t + T — y,m'')p{m'')drn: 
G(0; 0,t + T — y, rn'')p{rn")dm‘ 


'mo 

poo 


ft+T 


=n(<h(-) (g) A^(-, r))(f) + n J ^{y)N_{t + t - y)dy 
=n(<h(-) (g) A^(-, r))(f) + n f + r — x)N_{x)dx, 

Jo 

where we have used equations (155]) and fl39p and in addition, in third and last equality, 
we have made the change of variable y = t + r — x. It then follows that 

N(t, t) =n(<I>(-) ® N{a r))(f) + n{^{t + •) (g) N_{-)){t). (46) 
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Integrating both the members with respect to times, between zero and inhnity, we get 


^oo rt 


N{t, T)dt =n 


'0 ^0 


^oo nr 


+ n 


JQ JO 

r‘00 POO 


=n 


^{t — x)N{x, T)dxdt 

<h(t + r — x)N_{x)dxdt 
— x)N{x, T)dtdx 


'0 J X 


T roo 


n 



+ r — x)N _{x)dtdx 


'0 JO 

POO POO 

=n / dxN{x,T) / — x)dt 

Jo Jx 

PT POO 

+ n dxN_{x) / <h(t + r —x)(it 

Jo Jo 

POO PT POO 

=n / N{x,T)dx + n / dxN_{x) / ^{y)dy 

Jo Jo Jt-x 

POO PT 

=n / N{x,T)dx + n / N_{x)a{T — x)dx 

Jo Jo 

POO 

=n / A^(x,r)(ix + n(a(-) 0 A^_(-))(r), 

Jo 

where we have used equations (fT8|) and fl2n|) . We have then obtained equation (jUj). 

At this point, we are ready to proceed with the study of the approximation of 
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N_{t). Using equations fl26|) . flTTl) and fHTl) we get 

N_{t) =1 — '&[s(0; t, mo, m'),mo] + \I/[s(0; t, mo, m'),m] 


=1 

^1 


' mo 


J m 

poo 

- / p{m')p{m')s{0;t,mo,m')dm' 


' mo 


p{m')dm' — / p{m')p{m')s{0;t,mo,m')dm' 


p{m')p{m')s{0; t, mo, m')dm' + Q 


'mo 


where 


— / p{m')p{m')s{0;t,mo,m')dm' 

J m 

poo 

=Q + / p{m')p{m')dm'b{t) 

Jmo 

poo 

— / p{m')p{m'){^{-) D^{0; ■,mo,m')){t)dm' 

Jmo 

poo 

+ / p{m')p{m'){^{-) D_^.{0; ■,m,m')){t)dm' 

J mo 

poo 

— / p{m')p{m')dm'h{t) 

J m 

poo 

+ / p(m')p(m')(<l'(-) 0 Zi)+(0; •, mo, 

) 

p(m')p(m')($(-) 0 -D+(0; •, m, m')){t)dm', 

poo 

Q= p{m')dm'= 


' m 
poo 


(47) 


(48) 


As previously done, we study separately all the integrals in the final expression of fH7|) . 
Regarding the first three integrals, recalling again equations f[T5|) . flT6|) and flHOj) . we 
obtain the following results: 


p{m')p{m')dm'b{t) = nb{t)-. 


' mo 


p(m')p(m')(<h(-) 0 -D+(0; •, mo, m')){t)dm' 


' mo 


dm!p{m!)p{rn!) / dx $(t — x) G(0; 0, x, m^')p{m^'\m')dm'^ 


'mo 

rt 


' mo 


r*oo 


= / dx^{t — x) / c/m"G'(0; 0, x, m") / 


'mo 

poo 


'mo 

= / dx^{t — x) I G(0; 0, x, m")/A(m")(im"; 
v/ 0 J mo 
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and in the end, similarly, 


/ p{m')p{m'){^{-) -D+(0; •, m, m')){t)dm' = 

'mo 

pt poo 

= / dx^{t — x) / G{0-,0,x,m")lAi'm")dm" 


In the same way, for the last three integrals in W!\ we get the following expressions. 
The hrst one becomes 

poo poo 

/ p{m')p{m')dm'h{t) = b{t) / 

J m J m 

poo 

J m 

= (.-G-a)(m'-mp) 


-(/5 - a) 

^ ^Q-(P-a){m-mo) 

13 — a 


b{t) 


where we have used again equation ([TB]) and the ones (172]) and fl71|) in Appendix 
The second integral becomes 

poo 

/ p{m')p{m'){^{-) ® T*+(0; •, uiq, m')){t)dm' 

J m 

poo pt poo 

= / dm'p{m')p{m') / dx^{t — x) / G(0; 0, x, 

J m Jo J mo 

pt poo poo 

= / dx^{t — x) / dm"G{0-,0, x,m") / p{m')p{m')p{m''\m')dm' 

Jo J mo J m 

pt poo 

= I dx^{t — x) / G{0;0, x,m")lB{'m")dm", 

J mo 


>0 


where we have set 


= / p{m')p{m')p{m''\m')dm'. 


(49) 


In the end, the third one is similarly 

poo 

/ p(m')p(m')(<!)(•) (g) i7+(0; •, m, m')){t)dm! = 

J m 

pt poo 

= / dx^{t — x) / G(0; 0, X, 

Jo J m 

Substituting the obtained results in (1471) . we are then able to hnd the approximated 
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expression of the function N-{t): 

N_{t) + nb{t) 

pt poo 

— / dx^{t — x) / G{0-,0,x,m")lAi'm")dm" 


'mo 

poo 


+ / dx^{t — x) / G{0;0,x,m")lA{m")dm'‘ 

Jo J ni 

pt poo 

+ / dx^{t — x) / G(0; 0,x, 


'mo 

poo 


dx^{t — x) / G{0-,0,x,m")lB{m")dm'‘ 


=Q + nb{t) 


_ g-(/3-a)(m-mo) 


— / dx^{t — x) / ^(0; 0, X, 


'mo 

/*OD 


+ / dx^{t — x) / G{0;0,x,m")lA{m")dm" 


+ / —x) / G(0; 0, X, 


'mo 

/*oo 


(ix<f)(f —x) / G{0-,0,x,m")lB{m")dm'‘ 


( 60 ) 


At this point we can work on the approximation of the function A^_ (t) as we made 
for N{t, t). In particular, the hrst condition to impose is I Aim") = np{m"), justihed in 
the Observation [21 This is guaranteed by our proposal of triggered events’ magnitude 
conditional density function; 

p{m"\m') = p{m") [l + f{m', m")~\ 




1 + (7^(1 _ 2e-0-OK-™o)j ('x _ 2e-/5K'-™o)j 


for further details about this conditional density law see Appendix O On the other 
hand, in order to hnd the condition that the integral must satisfy, we substitute 
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in its expression, found in (jl9j), the above-mentioned proposal. It holds 


roo 

IsiiTi") = / p{m')p{rri')p{rn''\m')dm' 

J rn 


J m 

poo poo 

=p{m") / p{m')p{m')dm' + p{m'') / p{', 

J m J m 

poo 


m')p{m')f{m', m")dm' 


poo 

+ p(m") / p{m')p{m')Ci (l - (l - dm' 

J m 

[m—mo) 

poo 

_ / p(m')p(m')(l - 

J m. 


=p{m'')ne 


- p{m")Ci (l 


2g-(/3-a)(m'-mo 

J m 

poo 


°^^dm' 


[m-mo) 


J m 


-a)(m—mo) 


/? 

j^g-(/3-a)(m-mo) 


) +p(m")C'i(l -2e-^(’ 

0(m"-mo) j^g-2(/3-a)(i 

+ p(m")C'i(l - 2e-^(™"-”^° 
=p{m’')nH + ^(m")^! (^1 - 2^^ 


=p(m")n. 

+ p(m")Ci(l-2e-^('""-'"°) 


-2(/3 - 

-mo) I 

m 


a) (m—mo) 


/here the last equality follows observing 


/ riiL 

)j j^g-2(^-a)(m-mo) 

(5 


The expression obtained for integral iBim") depends then both on the Gutenberg- 
Richter function p{m") and on its square. This complicates a lot our study. However, 
observing that |1 —2e“^^™'"“™'°^| < 1, Gi < 1 and that the difference H — H^ is generally 
small for the typical values of parameters, we can approximate with a function 

of kind p{m")nL^ where the constant L is very close to H. As an example, we consider 
the following parameters data-set: 


13 = 2.0493, 
a = 1.832, 

Cl = 0 . 8 , 
mo = 1, 
m = 1.8. 
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x=Magnitude 


Figure 1: Plot of the functions f{x) = Hp{x) + Cip{x){l — 2£^)iJ(l — H) and g{x) = Lp{x), where 
p{x) = 13 exp{-/3(a: - mo)}. 


By means of a numerical algorithm minimizing the difference between lB{rri") and 
p{m”)nL, we get the constant L = 0.8032 that is, as we expected, close to the constant 
H = 0.8404. A graphical check is given in Fig. [1] 

We can then proceed with the substitution in equation (l50|) of the two conditions 

lAijn”) = np{m") 

iBifn") ~ nLpijn"), (53) 

for a certain constant L to be found numerically. 

Substituting again the symbol of approximation with equal sign for simplicity of 
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notation, we have 


N_{t) =Q + nb{t)[l - H] 


— n I dx^{t — x) I G{0;0,x,m")p{m")dm'' 


' mo 
poo 


+ n I dx^{t — x) I G{0-,0,x,m")p{m")dm" 


+ nL I dx^{t — x) I G{0;0,x,m'')p{m'')dm'' 


'mo 

poo 


— nL dx^{t — x) / G{0-,0,x,m")p{m'')dm' 

Jo J m 

=Q + nb{t) — nb{t)H + nLb{t) — nLb(t) 


— n I dx^{x) I G{0]0,t — x,m")p{m")dm" 


'mo 

poo 


+ n I dx^{x) I G{0;0,t — x,m")p{m")dm'' 


+ nL I dx^{x) I G{0]0,t — x,m")p{m")dni" 


' mo 
poo 


— nL dx^{x) / G{0-,0,t — x,m")p{m")dni‘ 

Jo J m 

=Q — nb{t)H + nLb{t) 


+ n[l — L] dx <h(a:) 


G(0; 0,t — X, ni")p{m")dm" 


mo 


(7(0; 0, t — X, m'')p{m")dm! 


J m 

=Q — nb(t) [H — L] + n\l — L~\ 


^{x)N _{t — x)dx, 


Jo 

where we have used equation ([39]). We then obtain 

N.{t) =Q - nb{t)[H - L] + Sm-) (54) 

where 

6 = n[l — L] . (55) 

Before proceeding with the research of the function N-{t), let’s do the following 
observation. 

Observation 5. If we integrate both the members of fl54|) with respect to time between 
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zero and infinity, recalling eqnations fl2n|) and fl2T]) we get 


r nt 


N_{t)dt =Qt — n[H — L] / h{t)dt + 6 

Jo 

=Qt — n[H — L] I 1 — a{t) 



0 Jo 


'0 


dt + 6 


— x)N _{x)dxdt 

— x)N _{x)dtdx 



=Qt — n[H — L] 
=Qt — n[H — L] 
=Qt — n[H — L] 


T — a{t)dt 


'0 J X 

+ 3 dxN_{x) / ^{y)dy 


T — A{t) 
T — A{t) 


+ S N_{x)[l — a{T — x)]dx 

Jo 

+ S [ N_{x)dx — 6{a{-) ® 


Hence, it is possible to find the convolution (a(-) ® A^_(-))(r) as a function of the 
integral N_{t)dt and vice versa. In particular, we hnd 


N_{t)dt 


1-6 


T 


Q-n[H-L] +n[H-L]A{T)-6{a{-)®N4-)){T)\ (56) 


and 


r 


{a{-)^N_{-)){T)=-[Q-n[H-L] 


n[H-L] , 

+ - -A{t) + 


1 - 


hJ 


N_{t)dt. (57) 


Using equation (I57)) . we can rewrite the function L(0; r) (equation only in terms 
of JJ" N_(t)dt: 


L(0;r) 


n 


1 — n 


Q-n[H-L\ n[H-L\^^^ ^-IT— ,,, 

-L-J + ^ ^ U (r) + ^ / N_{t)dt 


+ / N_{t)dt 


= / N_{t)dt 


n h — 1 
1 — n 6 


+ 1 


Q-n\H-L\ nAH-L] ,, , 


Recalling flM)) . let’s multiply equation fl58|) by the rate u relative to background events. 
We obtain 


a;L(0; r) 


where we have set 


5-1 


wA / N_{t)dt + T 


n u tiJ[H — L] 


+ OJ 


1 — n 6 5(1 — n) 


A = 


n 


1 — n 1 — 5 
and 

uj = uj\^Q — n[H — L]]. 


A{t), ( 59 ) 


( 60 ) 


( 61 ) 


If we could have an expression for fj N_{t)dt, using equation fl5^ in flM)) it would 
be possible to hnd the probability of having zero events in [0,r]. 
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In order to obtain an explicit expression of the function we use the standard 

technique of Laplace transform. More precisely, we integrate between zero and inhnity 
both the members of equation 0541) multiplied by the suitable exponential function. We 
get 


e-^^N_{t)dt= / e-^^Qdt - n[H - L] / e-%{t)dt 


+ 6 e-^\^-)®N_{-)){t)dt 

Jo 

= ^-n[H-L] [ e-**[l-a(t) 


dt 


+ 6 e-^^^{t)dt seC 


It then follows that 




1 — h e [ s 


+ re-‘am 


( 62 ) 


Since 


e ^^a{t)dt = 


^—st 


^—st 


^{x)dxdt 


1 — / ^{x)dx 


e-^^dt 


^—st 


dt 


^{x)dxdt 


Jo 

1 

s 


—st lOO 

lo 


^{x) e dtdx 


1 _i 
s s 
1 


'0 


1 - 


^{x)e ^^dx 
e~‘^^^{x)dx 


'0 


we obtain 


e-^^N_(t)dt 


s[l — 


Q-n[H-L\ / e-^^^{t)dt 


( 63 ) 


We notice that the above function of the complex variable s has one pole in s = 0. 
Furthermore, there are no poles with real part larger than zero. In fact, the other 
factor in the denominator is never zero for 3fJ(s) > 0, since 5 < 1 and 


e ^{t)dt < 


< 


'0 


< 

= 1 . 


e-^^^{t)dt 


\e-^^^t)\dt 

Mt)\dt 


0 
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Let’s then compute the Laplace transform of the function $(■) appearing in fl6^ . 


It is 


= 9c^ 
= 9c^ 


-i-e - 


X 


+ c 


-dx 


where 


= 9{scfe^^ / + sc)-^-^dx 

Jo 

poo 

= 9{scYe^'^ / e-yy-^-^dy 

J SC 

= 9{scYe^T{-9,sc), 

poo 

T{l,x) = / e~^t’'~^dt, |argx| < tt 


(64) 


is the incomplete gamma function (see USD. We notice that the integral in the third 
line of (164|) is the explicit form of the complex integral along the horizontal half-line 
starting from the complex point sc, written in the consecutive line. 

Substituting in equation fl65]) the solution (IMD . we obtain 


£[A^_(f)](s) = / e-^^N_{t)dt 

Jo 

1 


Q-n[H- L]9{scYe^^T{-9, sc) 


s[l — 69{sc)^e^^r{—9, sc)] 


(65) 


One possible way to hnd the function N-{t) consists of computing the inverse Laplace 
transform through the Bromwich integral; 


f*a+ZOO 


N(t) = - 

^ ^ 2m 


e 2.[N_{t)]{s)ds, 
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where a is any positive number. We then get 


^+oo 


N-{t) = 


27r 




—oo 
at /•+00 


=— / e'f‘£|iV_(i)l(<T + *{K 


' —OO 

^crt /•+00 


= ^ / ['2[iV-(^)](cr + ^0cos^^ + ^£[A^_(^)]((T + ^0sin^^]h^ 


— OO 

at ( r+oo r- 


e 


' —oo 

r' + OO p 


+v 

./ —oo 
( /•+00 ^ 


e 

+« 


3?(ii[A^_(t)](cr + i^)) + ?Q'(£[iV_(t)](cT + z^)) cos^td^ 

3?(ii[A^_(t)](cr + z^)) cos^t — Q'(£[A^_(t)](cr + Z(^)) sin^t 


' —OO 
^+oo p 


(£[iV_(t)](cr + z^) cos^t + 3fJ(£[A^_(t)]((7 + z^) sin^t d^ 


( 66 ) 


where 


^2.[N_{t)]{a + tO) =^yj N_{t)e-^^+^^'^*dtj=l N _{t) cos ^tdt 

and 




3 


{S.[N_{t)]{a + iO) =^[ I I =-y e-'"*iV_(t) sin^tdh 


Since A^_ (t) is a real function, we have 


r’+oo 


(il[A^_(t)]((T + z^)) cos^t + 3fJ(£[A^_(t)]((j + z.^)) sin,^t 


dC = 0; 


actually, this also follows easily from the odd-parity of the integrand. Moreover, ob¬ 
serving that the function 

W{^) = 3fJ(£[iV_(f)](cr -|- z^)) cos^t — Q'(£[iV_(f)](cr -I- z^)) sin^f 
is of even-parity, the last expression in fl66l) becomes: 


..crZ roo 


TT 


3fJ(il[iV_(f)]((j -|- z.^)) cos^t — Q'(il[A^_(f)]((T -|- z.^)) sin,^f 


d^. 


Finally, recalling that N-{t) = 0 for t < 0, one can verify that 


^+oo p 


3fJ(il[A^_(f)](cr -|- z^)) cos^t + Q^(il[A^_(f)]((j -|- sin,^f 


dc = o. 
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Then, 


3fJ(£[A^_(f)]((T + z^)) cos^td^ 

J Q-n[H- L]Ol(a + ,Oc]'’e<-*-«°r(-0, (a + »0c) \ 

Y (cr + ~ d9[{a + 0, (cr + *(^)c)] j 

(67) 

One may use a numerical approach to compute the last integral and then to hnd the 
required function N-{t). Once done this, we must integrat with respect to time. Then, 
the obtained result must be substituted in equation fl5^ . In conclusion, inter-event 
time density is obtained recalling that 

Frir) = 

A dr^ 

Because of the difficulty of both the computation of iV_(f) and the inter-event time 
density, we simplify the case we are analyzing by a further approximation. 

If we consider again r sufficiently small, we can set a(r — t)^ a{T) in (156]) . Substi¬ 
tuting once again the symbol of approximation with the one of equality for simplicity 
of notation, we get 


N.{t) = 


2e 


at 


TT Jo 

2^at roo 

TT Jo 


(1 — h) / N_{t)dt =T Q — n[H — L] + n[H — L]A{t) — 6 / a{T)N_{t)dt, 
Jo ^ ^ Jo 

from which we have 


N-{t)dt 


1 


1 — h -1- (5a(r) 


Q-n[H-L] +n[H-L]A{T) 


(69) 


The above-mentioned approximation makes sense since, if f G [0,r] and r in small, so 

is t. 
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We can then conveniently snbstitute the expression fl6^ of N-{t)dt in ([59]): 

1 


ujL(0;t) =^-— — uA 


1 — 5 + 5a(r) 


r 


Q-n[H-L] +n[H-L]A{T) 


1 T7 71 

+ Toj-- -h ~ 


h- 1 


61 — n h(l — n) 
uAt- 


6 1 — 5 + 6a{T) 

h-1 1 

H-— OjA 


1 — 5 + Sair) 


n[H - L\A{t) 


1 77 77 

+ TUJ-- -h — L\A{t) 


TU 


6 1 — n h(l — n) 
{6-1) A ^ n 


6 [l — 5 + ha(r) 1 — n 


+ 


un[H — L]A{t) 


{6-1) A ^ n 


TOJ + un[H — L]A{t) 
TU + un[H — L]A{t) 


1 — 5 + 6a{T) 1 — n 

5-1 

^ + 




n 


1 — 5 + 5a(r) (1 — n)6 

1 — n + na{T) 


(1 — n)[l — h + 5a(r)] 
where we have used equation fl6T]) and the last equality follows because 


(70) 


—A 

5 ^ 


+ 


n 


6-1 


1 — 5 + 5a(r) (1 — n)6 


n 6 

+ 


1 — n 5 — 1 


+ 


n 


+ 1 


1 — 5 + 5a(r) (1 — n)6 

n 

+ 


1 — 5 + 5a(r) (1 — n)6 


n 


n{6 — 1) 

(1 -n)5 

77,(5 — 1) + (1 — n)6 
(1 — ?7)5[1 — 5 + 5a(r)] (1 — n)6 

n{6 — 1) + (1 — n)6 + 77(1 — 5) + n6a{T) 
(1 — 77)5[1 — 5 + 5a(r)] 

1 — 77 + na{T) 


(1 — 77 ) [1 — 5 + 5a(r)] 


In this last calculus we have used equation (1^ . 

Now, using fl^S]) . we are able to hnd the probability density function Ft-{t) relative 
to inter-event time. From the calculations to obtain the second derivative of IP('r) = 
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e (equation fl9^ ). shown in Appendix iDl we get 


FJt 


_^g-a;L(0;T) , 


u + ujn[H — -h]a(r) 


+ ( rcj + ijjn[H — L\A{t) 
-ujn[H -L\^{t 


1 — n + na(r) 

(1 — ?7,)[1 — 5 + ha(r)] 


$(r)(h — n) 


(1 — n)[l — 5 + ha(r)]2 
1 — n + naij) 


+ 2 + ijjn[H — L]a{T) 
+ + ojn[H — L]A{t) 


(1 — n)[l — h + 5a(r)] 
<h(r)(h — n 


(1 — n)[l — (5 + (5a(T)]2 
(5 - n)<h(r)[(0 + 1)((5 - 1) + 5a{T){6 - 1)] 


=-exp<' - 


Tu + ujn[H — L]A{t) 
Zj + un[H — -h]a(r) 


(1 — n)(r + c)[l — 5 + ha(r)]3 

1 — n + na(r) 

(1 — ?7,)[1 — 5 + ha(r)] 

1 — n + na{T) 


+ ^ro; + un[H — L\A{t) 
-un[H - L\^{t) 


(1 — n)[l — h + ha(r)] 

$(r)(h — n) 

(1 — n)[l — (5 + ha(r)]^ 
1 — n + na(r) 


+ 2 + ijjn[H - L]a(r)^ — 

+ [too + ojn[H — L\A{t) 


(1 — n)[l — h + Sa{T)] 
<l>(r)(h — n) 


n)[l — (5 + ha(r)]2 

(5 - n)<h(r)[(0 + 1)((5 - 1) + Sa{T){9 - 1)] 
(1 — n)(r + c)[l — 5 + ha(r)]3 


(71) 


We notice that, differently from the case m = niQ, in this case m > tuq, in which 
not all the events are observable, our hypothesis of a conditional probability density 
function for the magnitudes of triggered events plays a role. In fact, this conditional 
distribution appears in expression flTTH through some of the constants, as for example 
L. 


5 Conclusions and future work 

We proposed and analyzed a new version of the ETAS model, where triggered event’s 
magnitudes are probabilistically dependent on triggering event’s ones. Thanks to the 
tool of the probability generating function and the Palm theory, we obtained a closed- 
form for the density of inter-event time for small values. More precisely, the closed-form 
is found by a linear approximation around zero of an exponential function. This is the 
only approximation needed for the case m = rriQ. In the general case m > mo, we have 
used an additional approximation to get the result in a closed-form. However, not using 
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the second approximation, we were able to find the studied density function in terms of 
the inverse Laplace transform of a suitable function. The results obtained show that our 
hypothesis of conditioning plays a role when we do not observe all the triggering events, 
as happens in practice. In fact, if we consider the particular case m = tuq, that is if we 
observe all the events that are able to trigger, the density obtained for the inter-event 
doesn’t depend on the new proposal for the distribution of triggered events’ magnitude. 
It depends only on the Omori-Utsu law and on its temporal integral. On the other 
hand, in the case m > ttiq, in which not all the triggering events are observable, we 
can see the influence of the above-mentioned proposal on some constants appearing 
in iv('r)- The inter-event time density then depends on the Omori-Utsu law, on its 
temporal integrals and on the new probability density function for triggered events’ 
magnitude. A possible problem to be analyzed for future work is the asymptotic 
study of the inter-event density for long times. Furthermore, it could be interesting 
to implement the new model and then to perform a simulation study. The statistical 
inference issue could be also of interest. 


Appendix 

A Phenomenological laws 

The Gutenberg-Richter law 

p(m') = 61 nl 0 - 
= &lnlO • 

_ (72) 

with (3 = 6 In 10, describes the empirical distribution of magnitudes. The value mo is 
the minimum magnitude an event must have to be able to trigger its own progeny. 
This law is adopted for events’ magnitude when we do not consider the respective 
characteristics of past seismicity. 

The Omori-Utsu law 

can be interpreted as the probability density function of random times at which first 
generation shocks independently occur, when the triggering event has occurred in t = 0. 
The productivity law 

p(m') = 

_ _„alnlO(m'—mo) 

— /Ac- 

= (74) 

where a = a In 10, represents the contribution of first generation shocks triggered by 
an event with magnitude equals to m'. 

For the previous three laws, the parameters (/9, 0, c, /«, a) are positive and typically 
estimated using maximum-likelihood algorithms. 
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B Probability generating function and some of its 
fundamental properties 

The probability generating function is a very useful tool to study random variables. 
It has several properties among which we can enunciate the following ones, which are 
very important for our analysis. 

Proposition 1. Let {Xi, X 2 , ■ ■ ■} be a sequence of independent and identically dis¬ 
tributed random variables, with generating function Gx- Let also N be a discrete 
positive random variable independent of Xi \/i, with generating function Gx- Then, 
given the random variable S defined as 

N 

i=l 

one has 

Gsi^) = Gn[Gx{z)]. (75) 

Proposition 2. Given two independent random variables X and Y, it holds 

Gx+y{z) = Gx{z)Gy{z). (76) 


C Conditional probability density function for trig¬ 
gered events’ magnitude 

The probability density function we propose for triggered events’ magnitude is not 
the Gutenberg-Richter law. Indeed, it is more realistic to suppose that the law of a 
triggered event’s magnitude depends on the one of the corresponding triggering event. 
We then consider a conditional probability density function with respect to triggering 
event’s magnitude. In what follows the latter magnitude will be indicated with m'. We 
search for a function p{-\m') such that: 

poo 

/ p{m''\m')dm'' = 1; (77) 

J mo 


p{m''\tr) = pfm”), 


(78) 


where 



p{rn')p{rn')p{m''\rn') , 

Imo t(^0 pi'm') dm' 


(79) 


indicates the probability for an event to have magnitude m", given the fact that it 
is not spontaneous. Recall that p{m') is the Gutenberg-Richter law (equation (172]) in 
Appendix 1X1) . p{m') is the productivity law (equation (17^ in Appendix lA]) and mo is 
the reference magnitude, that is the minimum value for an event to trigger other shocks. 
The above-cited condition fl78|) corresponds to the need of obtaining the Gutenberg- 
Richter law when we average over all the triggering event’s magnitudes, as previously 
explained in Observation [31 
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Observation 6. The right side of definition can be derived as below. Let’s consider 
the three following random variables: 

• M' ~ Gutenberg-Richter law is the triggering events’ magnitude; 

• N counts the number of shocks that a generic event triggers; it is such that 
p{N = k) = J^^p(N = k\m')p(m')dm'; 

• counts how many triggered shocks among the previous N have magnitude 
m'^ G [m", m" + h). 

Since 

( pm''+5 \ 

k, J p{x\m')dx\, 

where Binx{n,p) indicates the law of a binomial random variable X with parameters 
n and p, we have 

poo / pm"-\-5 \ 

p{^m'^ = q) = / / ^ Binx^„=q ( fe, / p{x\m')dx\p{N = k\m')p{m')dm'. 

Let’s consider the probability for a triggered event to have magnitude in the interval 
+ h). In order to find an expression of the just-mentioned probability, we 
suppose to have n realizations of the three variables We focus on the 

ratio between the sample means of the random variables and iV*, for i = 1,... ,n. 
On one side, this ratio converges by the law of large numbers to the above probability. 
On the other side, it also converges in probability to the ratio of the respective expected 
values, that are the quantities 


EJlo 9 XTo T,T=o (^k, p{x\m’)dx'^p{N = k\m')p{m’)dm' 

T.T=o k IZ = k\rn')p{m')dm' 

IZ EJlo qBinN^,,=q(^k, '^Z{x\m')dx^p{N = k\m')p{m')dm’ 

IZ ^kLo = k\m')p{m')dm' 

_ IZoT^kLo^ p{x\m')dxp{N = k\m')p{m')dm' 

fZ p{m')p{m')dm' 

_ IZo IZ" ^^P{x\m')dx J2 T=o kp{N = k\m')p{m')dm' 

IZ) Plkn')p{m')dm’ 

IZq IZ” p{x\'kri')dxp{rn')p{rn')drn' 

p{m')p(m')dm' ^ ^ 

Dividing by 6 and computing the limit for h tending to zero, we get precisely the 
probability for a triggered event to have magnitude m”: 

p{m''\m')p{m')p{m')dm' 
p{m"\tr) = ^ 
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where we recall that 

n = 

In equation (IHOj) we have used the fact that the productivity law p{m'), dehned in 
Appendix Rl (equation fl74p i. can be rewritten as 


p{m')p{m')dm! = 


(3k 


' mo 


/5 


p{m') = = k\m'). 

k=0 

In fact, it represents the mean number of hrst generation shocks triggered by an event 
with magnitude m'. 

We want to focus now on condition fl78p . which can be rewritten in the following 
way: 

POO 

/ p{rn') p{m')p{m''\rn')drri' = np{m"). (81) 

J mo 

Given the expressions of the functions p{m') and p{m'), let’s hnd a suitable law 
p{m"\m') such that flSip holds. Some qualitative behavior of this function has been ob¬ 
tained from the results of a statistical analysis of some Italian catalogues. The function 
we are looking for should be such that, when m' increases, the probability of having 
events with high magnitudes must increase and, at the same time, the one of events 
with lower magnitudes must decrease obviously. For further information, see |14) . 

In view of these considerations, the choice we adopt for p{m"\m') is 

p{m"\m') = p{m!') [l -|- f{m', m")], (82) 


with 

f{m', m") = —q{m') + s{m') -mo)'^ _ 

For now, the functions q{m') and s(m') are unknown. As we are going to see, they can 
be chosen in order to satisfy the following conditions: 

• p{m"\m') > 0, 

• p{m"\m') is such that condition (l8Tp is verihed, 

• p{m"\m') has the qualitative behavior previously introduced. 

Let’s proceed then with the search of the above-mentioned functions q{m') and s{m'). 
First of all, the conditional function p{m"\m') must be a density. We then have 


p{m'')[l + f{rn',m')~^dm''= 1 / p{rn'')f{m',m'')dm'' = 0. (84) 


'mo 


'mo 
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Substituting the expression chosen for f{m', m") we obtain 


0 = / p{m")[—q{'m') + s{m') — s{m')e 

Jmo 

pCX) 

= —q{m') + s{m') — s{m!) / 


I mo 


/ /\ / /X s(m') 

= —q{m) + s{m) --— 


= -q{m!) + 


s[m 


2 ’ 


from which 


s{m') = 2q{m'). 

Substituting the last result in equation (l8^ we get 

/(m', m") = q{m!) — 2q{m')e~^^'^ 


(85) 


Recalling equation fl82p . in order to have also p{m"\m') > 0, it should hold f{m', m") > 
— 1, that is 

g(m') - l) < 1. 

Since — 1| < 1, this is guaranteed if 


\q{m')\ < 1. 

Let’s impose at this point condition (I5T1) : 

poo 

/ p{m')p{m')p{m") [l + f{m\ m")\dm! = np{m") 

t 

r*00 

p{m!)p{m')f{m\m'')dm' — 0. 


( 86 ) 


' mo 


(87) 


'mo 


For f{m', m”) as in fl85p we get 

poo 

^mo 

poo 

= (l — 26“^*'™'”“™'°^) / p{m')p{m')q{m')dm'. 

Jmo 

A sufficient condition to solve the previous equation is obviously 

poo 

/ p{m') p{m')q{m')dm' = 0. 


( 88 ) 


' mo 


Before making a choice of the function q{m') such that condition (l88P holds, let’s make 
some qualitative comments. 
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The function q{m') should be such that the conditional law p{m"\m') has the qual¬ 
itative behavior before introduced. To this aim, we assume q{m') to be continuous and 
increasing. Furthermore we impose that q{m') is negative for m' < m at which it be¬ 
comes zero and positive elsewhere (see Fig. |5]). Let’s notice that when m' = m we have 
f{m',m") = 0 and the conditional probability density becomes the Gutenberg-Richter 
law. From equation fl85|) . it can be observed that, since the function 


is always increasing, the increasing or decreasing of f{m', •) will be determined by the 
sign of q{m'). In particular, for m' > m where q{m!) > 0, the function f{m',m") 
increases in m”. For smaller (higher) values of 


m" = — In 2 -|- mo, 
P 


where 1 — = 0, the function f{m\m") will have negative (positive) sign. 

Hence, the conditional density p{m"\m') = p(m")[l -|- /(m',m")] is below p{m") for 
magnitude values smaller than m" = ^ In 2 -|- mo and above it for higher magnitude 
values. Vice versa, for m' < m where g(m') < 0, the function f{m',m") is decreasing 
in m". With analogous reasonings we come to the conclusion that p(m"|m') is above 
the Gutenberg-Richter law for values smaller than m" = ^ ln2 -|- mo, below for higher 
values. 

Observation 7. For triggered events’ magnitude values such that m" > ^ln2 -|- mo, 
the factor 1 — is positive. Hence, hxed m" larger than the just-mentioned 

value, since q{m') is increasing, the function f{m', m") is increasing with respect to the 
first argument, too. This means that when triggering event’s magnitude increases, so 
does the probability density function at a point m" > ^ln2 -|- mo. On the contrary, 
the value of the probability density decreases at a value m" < ^ In 2 -|- mo. We are then 
in qualitative agreement with the graphical results obtained. 

To sum up, the function q[m') must be: 


• less than one in absolute value. 


• increasing, 

• negative for magnitudes less than a certain value m, positive otherwise, 

• such that condition fl88|) is verified. 

The choice we make for g(m') is the following: 

q{m') = -Gi + Ga(1 - g-O-OK-^no)) ^ (gg) 

where the constants Ci and C 2 are obtained imposing the above-cited conditions for 
q{-). We specify that in practice /3 > a: it is the condition which ensures that the 
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Figure 2: Plot of the function (/(a;) = (71(1 —2exp{—0.4648(a:—1.8)}). The value at which the function 
becomes zero is m* = 3.2913. 


process doesn’t explode. Firstly, let’s consider equation fl 88 |) . Recalling the expressions 
of p{m') and p{m^) (respectively equations fl72l) and fTMl) in Appendix!^, we get 


0 = (3k 


„-(/3-a)(m'-mo) 


'mo 


-Cl + C'2(l - 


dm' 


= (3k 


{-Cl + C2) / -C2 


' mo 


' mo 


= (3k 


-C 1 + C 2 


C 2 


(3 — a 2{(3 — a) 


(3 — a 

from which follows that 

Hence 


^2 = 2^1. 


q{m’) = Ci- = Ci (l - . (90) 

In order to have also \q{m')\ < 1, since |1 — < 1 , it is enough to choose 

1 > (7i > 0. One can easily verify that the function q{-) is increasing since Ci > 0 and 
(3 > a (Fig. [2]). Furthermore, if Ci 7 ^ 0, the value for which q{m') = 0 is 


m = m = 


(3 — a 


In 2 + mo. 


Since q{m') increases, it will be positive when m' is larger than this value and is negative 
elsewhere. 

At this point, substituting the expression flQUll of q{m') in fl85|) . we get 

f{m', m") = Ci{l- (l - . (91) 

In Fig. [3](a) we can see the plot of the function f{m', ■) obtained with the values of 
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(a) 



(b) 

Figure 3: Plots of the functions f{m',x) = 0.8(1 — 2 exp{—0.4648 (to' — 1.8)})(1 — 2 exp{—1.9648(x — 
1.8)}) and f{x,m") = 0.8(1 — 2exp{—0.4648(a; — 1.8)})(1 — 2exp{—1.9648(m" — 1.8)}), 
respectively in (a) and (b). In the first subplot the value that varies is the triggering events’ 
magnitude. On the other hand, in the second one the varying value is the triggered events’ 
magnitude. 


parameters such that we get + tuq = 3.2913. This is the value at which q{m') be¬ 
comes zero. The values of parameters correspond to a real situation because they were 
estimated based on a real data set. As expected, decreases (increases) with 

respect to m" when m' is smaller (higher) than the just-mentioned value. Furthermore, 
it holds f{m', m”) = 0 for m” = ^ + = 2.1527. In Fig. I2^b) one can observe that, 

for any hxed value of m" > 2.1527, the function is increasing with respect 

to the hrst variable. 

In conclusion, substituting equations and fl72|) (Appendix]^ in fl82l) . we obtain 
the probability density function relative to triggered event’s magnitude, conditioned 
on the magnitude m' of its own triggering event: 






1 + ( 71(1 - (1 - , 


(92) 
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for all 1 > Cl > 0. 

We study now the behavior of p{m''\m') by looking at its derivative. Recalling 
equations (15^ and (153]) we have 


df{m\ m!' 


dm" 


dm" 

= —/3p{m")[l + f{m', m")] + p{m") 

= /3p{m") [-1 - f{m',m”) + 

= /3p{m") [-1 - g(m') + + 2g(m')e-^(”^"-™° 

= f3p{m") [—1 — q{m') + . 


It follows that if q{m') < 0, that is m' < ■g^ln2 + mo, since \q{m')\ < 1, the density 
function p{m"\m') is always decreasing in m". If instead q{m') > 0, that is m' > 
■gWiii2 + mo, the above-cited density increases in m" till a certain maximum reached 
in 

= 1 -|- qim') m" = \ In ^ -|- mo 

^ ^ ^ ^ 13 l + q{m') 

and after that it decreases. This kind of trend is shown in Fig. 0] This behavior 
agrees with the results of the analysis we made in [T3]. In Fig. 0] we can also verify the 
behavior described in Observation [T] 
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x=Magnitudo 


Figure 4: Plot of the probability density function relative to triggered events’ magnitude when trig¬ 
gering events’ magnitude varies and Ci = 0.8. 


D Time derivatives of the probability P(t) of having 
zero events in [0, r] 

In this appendix we include the calculations we made to obtain the second time deriva¬ 
tive of the probability P(t) = in the case m > uiq. This derivative is indeed 

useful in order to get the inter-event time density TV(r) by means of equation fl2^ . 
Using equations (12T|) and (TTn]) we have 




u + un[H — L]a(r) 


+ 


Tu + ujn[H — L]A{t) 


1 — n + nair) 

(1 — ?7,)[1 — d -b ha(r)] 
dr 1 — n + nair) 


dr i{l — n)[l — 6 + ha(r)] 


It follows that 


dP(r) 

dr 




+ 


(jj -b ojn[H — T]a(r) 
d 


TOO + un[H — L\A{t) 


dr 


1 — n + nair) 

(1 — n)[l — d -b Sa{T)] 

1 — n + na{T) 


(1 — n)[l — h -b da{T)] 
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The second derivative of P(t) with respect to r is then 


dr^ dr 


+ 


Tu + un[H — L]A{t) 


u + ujn[H — T]a(r) 
^ d 


1 — n + naij) 


1 — n + na(r) 


-e-^L(o,r)^ -a;n[i/-L]$(r) 


dr \_{1 — n)[l — 5 + J 
1 — n + na(r) 


[(1 - n)[l - h + (5a(r)] J 


+ 

+ 

+ 


u + ujn[H — -h]a(r) 
u + ujn[H — T]a(r) 

TU) + u)n[H — L\A{t) 


d 

1 — n + naij) 

dT 

(1 — n)[l — h + 5a(r)] 

■ d 

1 — n + na(r) 

dT 

(1 — n)[l — h + 5a(r)] 


1 

1 — n + nair) 

1 

dT‘^ 

(1 — n)[l — d + ha(r)] 

J 


_g-wL(0;r) ^ 


UJ + ujn[H - L]a(r)^ ^ 


1 — n + na(r) \ 


n)[l — h + 5a(r)] J 


+ { TU + un[H — L]A{t) 
-un[H - L]$(r 


$(r)(5 — n) 


(1 — n)[l — h + (5a(r)]2 
1 — n + nair) \ 


+ 2 + un[H - L]a(r)^ — 

+ (tu + un[H — L]A{t) 


(1 - n)[l - h + (5a(r)]/ 
$(r)(h — n) 


n)[l — (5 + ha(r)]^ 

(d - n)<h(r)[(0 + 1)((5 - 1) + ha(r)(6' - 1)] 
(1 — n)(r + c)[l — (5 + ha(r)]^ 


(93) 


The previous formula has been obtained using the following results: 


daij) 

dT 


A. 

dT 


1 — / ^{x)dx 


d 

1 — n + nair) 

d 

5 — n + n[l — d + ha(r)] 

dT 

(1 — n)[l — h + da(r)] 

dT 

(5(1 — ?7.)[1 — (5 + (5a(r)] 


5 — n d 1 

h(l — n) dT 1 — 5 + ha(r) 
$(r)((5 — n) 

(1 — n) [1 — h + 5a(r)]2 ’ 
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in the end 



1 — n + na^r) 

S — n d 

<l)(r) 

dr'^ 

(1 — n)[l — h + (5a(r)] 

1 — ndr 

[1 — 5 + ha(r)]2 


6 — n 


(1 — n)[l — h + Sa{T)Y 


—d>(r) _^ [1 — (5 + ha(r)]^ 


r + c 


+ 2<h(r)[l — h + 5a(r)](5<l>(r) 


(h — n)$(T) 


0 + 1 

(1 — n)[l — h + (5a(r)]^ ] r + c 
(5 — n)<l)(r) 


(1 — n)(r + c)[l — 5 + ha(T)]3 
{6 - n)<l)(r)[(0 + 1)((5 - 1) + 5a(r)(0 - 1)] 
(1 — n)(r + c)[l — h + 5a(r)]3 


[—1 + 5 — Sa{T)] + 2h$(r) 

(0 + l)(h — 1) — S9a{T) — 6a{T) + 250a(r) 


The last computation has been got since, recalling equation (120]) and the one fITO]) in 
Appendix]^ we have 


and 


(i<h(r) 

dr 


(jj'T 

-0c^0+l)(r + c)-"-2 

0 + 1 


-d)(r) 


r + c 


$(r) 


0a(r) 
r + c 


45 




















Notation 


a, a 

exponent of the productivity law (equation (1741). Ap¬ 
pendix!^; it holds a = a In 10. 

a{t) 

auxiliary function defined in fl20|) as the integral of the Oniori- 
Utsu law (equation d73|). Appendix!^ between t and infinity. 

A{t) 

auxiliary function defined in 0211) as the integral of a{x) be¬ 
tween zero and t. 

A 6 

exponent of the Gutenberg-Richter law (equation fl72|). Ap¬ 
pendix 1^; b is the so-called b-value and holds that (3 = 
61n 10. 

b{t) 

auxiliary function defined in (jH]) as the integral of the Oniori- 
Utsu law (equation (175]). Appendix!^ between zero and t. 
It holds bit) = 1 — ait). 

c 

parameter of the Omori-Utsu law (equation flTHj). Ap¬ 
pendix!^. 

Cl,C2 

constants entering in the definition of the function g(m'); it 
holds C 2 = 2(71. 

D{z] t, T, rriQ^rn) 

auxiliary function defined in fl2^ entering in the definition of 
yiz]t,T, mo,m'). 

t, fh,m) 

auxiliary function defined in fl30|) entering in the definition of 
yiz; t, r, mg, m') and siz; t, niQ, m'). 

5 

A 

f{m', m”) 

parameter defined in 0551). 
parameter defined in flOOl). 

auxiliary function defined in (Appendix !C|) useful to cal¬ 

culate the conditional probability density p(m" m'). 

f{z, m', m) 

function defined in (!2]) useful to calculate both the probability 
generating function of the total number of events in [0, r] and 
the one relative to the number of event triggered by a generic 
shock in the same time interval. The multiplication by this 
function allows us to add the triggering event if its magnitude 


is larger than the observability threshold m. 

Omori-Utsu law (equation (175]). Appendix [A!) : it can be in¬ 
terpreted as the probability generating function of random 
times at which first generation shocks independently occur, 
triggered by an event occurred in f = 0. 

FAr) 

G{z] t, T, m') 

inter-event time density. 

probability generating function of the number of observable 
events triggered in [0, r] by a spontaneous shock of magnitude 
m' occurred in —t, with f > 0. 

H 

I Aim") 

parameter defined in fl52l). 

magnitude integral from the reference cutoff mg to infin¬ 
ity, defined in flTSl). of the product between the Gutenberg- 
Richter law fl72]) (Appendix !^, the productivity law (1741) 
(Appendix !^ and the conditional probability density func¬ 
tion (Appendix !C]) of triggered events’ magnitude. We 

impose that I Aim") = np(m"). 
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K 

L 

L{z]t) 


X 

m 

mo 

n 

N.{t) 

u 

u 

Vt{z-, t) 
p{m) 


P(t) 

p{m"\m') 

p{m''\tr) 

PGF 
'^{h, m) 

Q 


q{m'), s{m') 


p{m) 


s{z]t,mQ,m'), 
y{z-,t,T, mo,m') 

e>o 


magnitude integral from the completeness value m to infin¬ 
ity, defined in of the product between the Gutenberg- 
Richter law fl72|) (Appendix |^, the productivity law flTH) 
(Appendix and the conditional probability density func¬ 
tion (1^ (Appendix [C]) of triggered events’ magnitude. We 
impose that /^(m") ~ nLp{m"). 

multiplicative parameter of the productivity law (equa¬ 
tion (jUj), Appendix [A|) . 

parameter entering in the approximate condition that the 
integral /^(m") must satisfy; it has to be found numerically, 
auxiliary function useful to compute the probability gener¬ 
ating function VL{z] t) and then the probability to have zero 
events in [0, r]. 

rate of the whole process of ’’observable” events, 
threshold value of completeness magnitude, 
threshold value of reference magnitude, 
average branching ratio dehned alternatively in fll6|) . 
auxiliary functions useful to compute L(0;r); it holds 
N{t,T) = Nm=moit,T). 

rate corresponding to the ’’observable” events of the sponta¬ 
neous component of the process, 
parameter dehned in (I6T1) . 

probability generating function of the total number N{t) of 
observable events in [0,r]. 

Gutenberg-Richter law (equation (I72|) . Appendix]^; it is the 
probability density function of events’ magnitude when we do 
not consider past seismicity, 
probability to have zero events in [0,t]. 
conditional probability density function assumed for trig¬ 
gered events’ magnitude; it is dehned in fAppendix iGl) . 

conditional probability for an event to have magnitude m" 
given the fact that it is triggered; it is dehned in fl79|) . 
Probability Generating Function. 

auxiliary function dehned in fl26p useful to compute L{z]t). 
parameter dehned in fHSll and calculated as the integral of 
the Gutenberg-Richter law between the completeness value 
m and inhnity. 

auxiliary functions entering in the dehnition of f{m', m") 
(equation flM]) h the function g(m') is dehned in (1^ and 
holds that s{m') = 2q{m'). 

productivity law (equation ([71]), Appendix]^; it is the con¬ 
tribution of hrst generation shocks triggered by an event with 
magnitude m'. 

auxiliary functions used to dehne L{z]r). 

exponent of the Omori-Utsu law (equation fl73|) . Ap¬ 
pendix]^. 
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